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We study how an oscillating mirror affects the electromagnetically induced transparency (EIT) of 
an atomic ensemble, which is confined in a gas cell placed inside a micro-cavity with an oscillating 
mirror in one end. The oscillating mirror is modeled as a quantum mechanical harmonic oscillator. 
The cavity field acts as a probe light of the EIT system and also produces a light pressure on the 
oscillating mirror. The back-action from the mirror to the cavity field results in several (from one to 
five) steady-states for this atom-assisted optomechanical cavity, producing a complex structure in 
its EIT. We calculate the susceptibility with respect to the few (from one to three) stable solutions 
found here for the equilibrium positions of the oscillating mirror. We find that the EIT of the atomic 
ensemble can be significantly changed by the oscillating mirror, and also that the various steady 
states of the mirror have different effects on the EIT. 

PACS numbers: 42.50.Tx, 03.67.Bg, 32.80.Qk 



I. INTRODUCTION 

Fast developments are now occurring in studies at the 
interface between different kinds of physical systems. Ex- 
amples of these include: couplings between light and 
nano/micro-mechanical systems (e.g., [J [2, H, 0, S S 
0)11, H, called optomechanical systems; interactions 
between superconducting artificial atoms (e.g., charge 
or flux qubit) and transmission line resonators (e.g., 
[H [H 111, Q E, 13), and so on. These studies are 
partly motivated by possible physical implementations 
of quantum information processing, to explore potential 
future devices, and to study the interesting physics in 
these hybrid structures. 

In optomechanical systems (e.g., P, [3, H, 0, H, 0, 0, H, 
fioj). the radiation pressure acts on the oscillating mir- 
ror and induces the interaction between the mechanical 
system and the optical field. The back-action from the 
mirror to the cavity field [Tt'I can result in several (about 
5) steady-state solutions for the equilibrium position of 
the mirror. It has been proved [l| that this system ex- 
periences bistability in some parameter region. It has 
also been shown that light pressure can be used to real- 
ize entanglement between the cavity field and a micro- 
scopic object (e.g., a moveable mirror) H 0, H, and 
can also cool down the mirror [1, H, 0, [8|. Furthermore, 
Refs. [H, [l^ studied optomechanical cavities contain- 
ing a two-level atomic ensemble. The atomic ensemble 
effectively enhances the radiation pressure of the cavity 
field on the oscillating mirror, producing a cavity-atom- 
mirror entanglement [18, 19]. In this paper we consider 
this atom-optical system with a more complex atomic 
ensemble, which can enable quantum interference, e.g., 
electromagnetically induced transparency. 

Electromagnetically induced transparency (EIT) is a 
remarkable quantum interference phenomenon, which en- 



ables the active control of light propagation in a coher- 
ent medium [13, [U, [H, HI]. Usually, the basic popu- 
lation transfer configuration for the atoms in EIT is of 
A-type, where the two transitions from a common up- 
per energy level to two lower energy levels are induced 
by two different optical fields (e^.jdassical control field 
and probe quantized field) 0, [2^, , respectively. One 
is a strong field, and the other is a weak one. The 
stronger field can effectively modify the susceptibility of 
the medium so that the weak one (as a probe signal) 
can pass through the medium transparently at the two- 
photon-resonance [13, m, m, HH. Recently, it has 
been shown that this effect can also work well at the sin- 
gle photon level for the probe light [s^ , and thus the weak 
field must be treated quantum mechanically. In the quan- 
tum approach, a dark state with dressed photons can be 
invoked to store quantum information of photons on the 
atomic ensemble as quantum memory (ssl . Is^ . These 
quantum manipulations at the single photon level re- 
quire frequency-matching with extremely high-precision 
for one- or two- photon resonances. When the quantum 
field is provided by a micro-cavity with a one-end oscillat- 
ing mirror, the oscillation of the mirror might affect such 
a precise frequency-matching condition and thus affect 
the EIT. 

Motivated by these works, here we study how the 
oscillating mirror changes the properties of the EIT. 
We will study the atom-assisted optomechanical micro- 
cavity, through which we explore the possibility to inter- 
face other systems, via some physical mechanism, such 
as EIT. Here, the atomic ensemble for each atom with 
A-type transitions is placed inside a cavity with a one- 
end oscillating mirror. Due to the mirror's oscillation, 
the susceptibility of this atomic medium displays a multi- 
stability corresponding to the multi-equilibrium positions 
of the mirror. Another prediction of our study is that 
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FIG. 1: (Color online) Schematic diagram for the atom- 
assisted optomechanical system considered here. There are 
three main components: (i) an optical cavity with one fixed 
mirror; (ii) another oscillating mirror, which is modeled as a 
quantum-mechanical harmonic oscillator (denoted by a black 
spring); and (iii) an ensemble of identical three-level atoms, 
which are confined inside the cavity. Each atom is assumed 
to have A-type transitions. Here, Ap — LUa — l^o + i^o (a^) /' 
and Ac = u>a — u)c- 



the mirror's oscillation significantly alters the properties 
of both the real and imaginary parts of the susceptibil- 
ity. We also investigate in detail how the different steady 
states of the mirror influence the dispersion relations and 
absorption properties of the light. 

The paper is organized as follows: in Sec. II, we intro- 
duce the model for the optomechanical system with EIT; 
in Sec. Ill, we present the Heisenberg-Langevin equations 
for this system and obtain several (from one to five, de- 
pending on the system parameters) steady-state solutions 
for the position of the oscillating mirror; in Sec. IV, we 
study the EIT with susceptibilities for different (from one 
to two) equilibrium positions of the mirror; in Sec. V, we 
summarize our results. 



II. ATOM-ASSISTED OPTOMECHANICAL 
SYSTEM 



a classical light field with frequency v induces a transi- 
tion between \a)^ and |c)j, which is often used as a control 
field. The quantized cavity field, with frequency wq when 
the oscillating mirror is fixed, induces the transition be- 
tween |a)j and \b)^. We assume that this cavity field acts 
as a probe field. Here, the transition between |6) . and 
|c)^ is assumed to be forbidden. The oscillating mirror 
is modeled as a quantum-mechanical harmonic oscillator 
with frequency ujm and the mass M . This harmonic os- 
cillator can also be considered as a spring with an elastic 
coefficient Mu}\j. 

Based on the above considerations, and using ft = 1, 
the Hamiltonian of the total system 



H — He + Hm + Ha + Hm-c + Ha-l 



(1) 



has five terms corresponding to the cavity {He), the os- 
cillating mirror {Hm) of mass M, the atom gas {Ha), the 
mirror-cavity {Hm-c), and the atom- light term {Ha-l)- 
Explicitly, these are described below: (i) the free Hamil- 
tonian 



H, 



c = ijjQO^a 



(2) 



of the single-mode cavity field with the annihilation and 
creation operators a and a^; this term (and wq) refers to 
a cavity with two fixed mirrors; (ii) the free Hamiltonian 



(3) 



of the oscillating mirror, where p is the momentum of the 
oscillating mirror with a small displacement x; (iii) the 
free Hamiltonian 



N 



Ha = y^Xuja(T 



(i) 
aa 



(4) 



of the N three-level A-atoms with the operators aaa = 
\a)ii{a\ and a = a, c; here, uja (^c) is the energy level 
spacing between |a)- and \b)- (|c) ■ and |5)-) for the ith 
atom, and we have assumed the ground state \b)i as an 
energy reference point; (iv) the interaction Hamiltonian 



H. 



M-C 



t 



(5) 



between the cavity field and the oscillating mirror 
presents a radiation pressure on the mirror due to the 
small change x of the cavity length when the mirror os- 
cillates where I is the cavity length when the mirror 
is at its equilibrium position; (v) the interaction Hamil- 
tonian 



As shown in Fig. [1] we consider an ensemble of N 
identical three-level atoms, which are confined inside an 
optical cavity with a one-end oscillating mirror. The ex- 
cited, meta-stable, and ground states of the ith atom are 
denoted by |a)j, |c)j, and \b)^. Each atom is assumed 
to have A-type transitions. That is, for the ith atom. 



Ha- 



N 

E 



(6) 



between the three-level atoms and the classical as well as 
the quantized fields. In Eq. ([61), is the Rabi frequency 
associated with the coupling between the classical field 
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and the three-level atoms. The frequency v is assumed 
to satisfy the condition v — w^ — Wc^^c- Here, is the 
detuning between the frequency of the classical control 
field and the transition frequency from the energy level 
\a)i to the energy level |c)i for the zth atom. 
The parameter 



5 = -/i^cjo/SFeo 



(7) 



in ([S]) describes the coupling between the quantized cav- 
ity field and the three-level atoms, where /i is the electric- 
dipole transition matrix element between levels |a) and 
1 6), V describes the volume of the cavity, and eg is the 
permittivity of the vacuum. We note that the effect of the 
oscillating mirror on the coupling between the atoms and 
the quantized cavity field [35| has been neglected when 
the Hamiltonian in Eq. ([1]) was derived, because we do 
not consider the strong coupling between the quantized 
field and the atoms. 

In the derivation of the Hamiltonian in Eq. Ilj , we have 
assumed that the linear size of the atomic ensemble is 
much smaller than the wavelengths of the light fields. In 
this case, the couplings between the atoms and the light 
fields are homogeneous, and we can define the collective 
operators of the atomic ensemble as in Refs. ^25,, jiQj 

AT N 

i—l i—1 

N N 

i=l i=l 

Together with Eq. ([8|), the interaction Hamiltonian Ha-c 
in Eq. ^ can be rewritten as 



H 



A-L 



= VI e-""* T+ + g\/7V a -|- H.c 



(9) 



We assume that the number N of atoms is large enough 
so that the collective operators in Eg. (El) satisfy the com- 
munication relations as in Refs. [l^ [25r[30j 



and 



[C\S]^Q, [A,S]=A, (10a) 



[r_,Ct]=0, [T+,C^]=A\ (10b) 



when most atoms are in their ground states. Equa- 
tions (jlOap and (jlObp present a dynamical symmetry in 
our system described by the semidirect-product algebra 
containing the algebra SU(2) with its generators T±. It is 
easy to prove, as in Ref. [19], that the collective operators 
and the communication relations can also be given in a 
similar way as in Eqs. ([6]), (|10a[) and (jlObp for the case 
when the couplings between different atoms and the light 
fields are inhomogeneous. Therefore, our study here can 
be generalized to the case of inhomogeneous couplings. 



III. HEISENBERG-LANGEVIN EQUATIONS 
AND MULTI-STABILITY 



A. Steady-state positions of the moveable mirror: 
analytical results 



Using the commutation relations in Eqs. (I10ap - (ll0bp . 
we can write the Heisenberg-Langevin equations of mo- 
tion as 

P 



dtx 



(11a) 



dtp = ~^P + y - MujIj X - V7l7ein {t) , (lib) 



dta = - ^ i.^^j)"'^ igVN A + y/%ain {t) , 

(11c) 

dtA = --f^A-iuJaA-ine-^^^C-igVNa+fiit), (lid) 

dtC = -72C - iojcC - in e""' A + f2{t). (lie) 

Here, dt denotes a time derivative. We have phenomeno- 
logically introduced the damping 7^^ of the oscillating 
mirror, the decay rates 7q for the cavity field and (72) 
for the decay from \a) to |6) (|c) to \b)), respectively. We 
also assume that the quantum fluctuations of the cavity 
field, mirror, and the atoms satisfy the conditions 



and 



(ei„ (i)) = (A it)) - (A (t)) = 



(Cin (t)) = Oin (t) . 



(12) 



(13) 



Here, ain (t) can be understood as an input driving field. 
That is, ain (t) can be rewritten as 

ain (i) = "in (t) + Sain (t) , 

where the quantum fluctuation of the input held Sai^ (t) 
satisfles (<5ain (t)) = 0. 

To obtain the steady-state solutions, let us flrst remove 
the fast varying factors by the following rotating trans- 
formations 

a = a exp(— iciJit), (14a) 
A = Acxp{-iujLt), (14b) 
C = Cexp[i{Ap- Ac- ujc)t], (14c) 

and 

Oin {t) = Oin (t) exp{~iujLt) (14d) 
= [ain{t) + Sain{t)]exp{-iujLt), (14e) 

where the detuning Ap between the transition frequency 
uja of the atom, between the energy levels \a) and \b), and 
the effective frequency of the cavity field, is given by 



Ap = - COL, 



(15) 
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(16) 



with the effective cavity frequency 



ujl=ujo J- {x} 



Here, (x) denotes the mean value of x. Equation pi 
shows that the effective frequency of the cavity can be 
changed by the oscillating mirror. When the osciUating 
mirror has zero displacement, then lol equals to uq. 

In the rotating frame, the Heisenberg-Langevin equa- 
tions in Eqs. (|llb[) - (|llep become 



dtp = -^P + y - Miult X - ^ein (t) , (17a) 



dtA^ -(7i + iAp)A - iflC - igV N & + fi{t), (17b) 



dtC = - [72 + * (Ap - A,)] C-tnA + hit), (17c) 



:^-t^{{x)~x) a~igVNA + ^ai,,it), 



dtCL = - 

where the fluctuation forces are 



/i(t) = /i(i)exp 



iujQ { 1 + — I t 



and 



f2{t) = h{t) exp [~i (Ap - A, ~ ujc) t] . 



(17d) 
(18a) 

(18b) 



We are interested in the steady-state regime. Let us first 
assume that all the time derivatives of the mean values 
for the operators in Eqs. (|17ap - (jl7d[) are equal to zero; 
then we obtain the steady-state equations 



'^{ii^)^{d)^~MLol,{x) 



0, 



(19a) 



^(a),-zgViV(A) +V7^ai„ = 0, (19b) 



(7i + iAp^s) (i)^ - in (c)^ - igVN (5), = 0, (19c) 



- [72 -t- z (Ap,, - Ae)] {^Cj^ - in ^Aj^ - 0. (19d) 

Here, (O) ^ (O represents the operator in the above 
steady-state equations) is the mean value of the opera- 
tor O under the steady state. The parameter Ap,, is the 
detuning described in Eq. (|15p when the system reaches 
steady-state. In Eqs. (|19ap and (jl9bp . the correlations 
(a^a)s and {xa)s have been approximately replaced by 
{a^)s{a.)s and (a;), (a),, respectively. These approxima- 
tions indicate that the correlations between the fluctua- 
tions near the steady states are very small compared to 



the corresponding mean values in the steady states. This 
can be quantitatively shown as (36| 



« 1, 



iiSx) (Sd))^ 



< 1. 



(20) 



From Eqs. (Il9b[) - (ll9d|) . when the system reaches the 
steady-state, the mean values (^4), and (a), can be easily 
derived as 



and 



A 



2a;, 



1 - 



-igy^NjQ Qin^(Ap,,) 

G(Ap,,) n{Ap^s) + ^n^' 
g^N n{Ap,s) 



with the functions 

il(Ap,,) 72 + « (Ap, 
and 



G(Ap,,) n{Ap^s) + ^n^ 

Ac) 



G(Ap,,) =^2^+^70 (7i+iAp,,). 



(21) 



(22) 



(23) 



(24) 



Using Eq. (I19a|) and Eq. ([22|) . we can self -consistently 
derive a nonlinear implicit equation for Ap,, as below 



1 - 



g^Nn{Ap^, 



G(Ap,,) l](Ap,,) + W 



7^ (Ap. - Ao) , 



(25) 

where the elastic constant (which has units of energy) of 
the spring attached to the moveable mirror is 



2 ;2 



and 



Aq = COq - LJo- 



(26) 



(27) 



for the detuning Aq between uja (the highest frequency 
of the A-atom) and ujq (the cavity frequency in Eq. ([2]) 
when the two mirrors are fixed). Equation (j25p is a 
fifth-power implicit equation for Ap,,. Thus the system 
may have several solutions (i.e., multi-stability in some 
parameter regions, corresponding to several steady-state 
positions of the mirror). These mirror positions are de- 
termined by 



(28) 



B. Steady-state positions of the moveable mirror: 
numerical results 



Let us first analyze the atom-cavity detuning Ap,, and 
the equilibrium positions (x) of the mirror. In principle. 



5 




FIG. 2: (Color online) Schematic diagram for the functions of 
the left Yl(Ap,s) and the right Yr(Ap,3) hand sides of Eq. pS)) 
versus Ap,s in (a) for a large value of Ao and (b) for Aq in 
the region of (a) indicated by the red dashed box. In (b), the 
brown triangle, red crosses, green dot and blue squares de- 
note four different steady state solutions of the self-consistent 
equation (|25p . The brown triangle (far right) denotes the solu- 
tion Ap^l, the four blue squares (top left) denote the solution 
Ap'^s , the red dots on the curve denotes the solution A^l , and 
the green dot denotes the solution A^^l. The solution Ap|l is 
not shown in the regions c and d. Here, A[1^ (i — ah, be, cd) 
denote the threshold values of Ao that separate regions with 
different number of solutions. Namely, A[J^''' is the boundary 
point of the regions a and b, Ay^'^' is the boundary point of 
the regions 6 and c, and Ajjj'*' is the boundary point of the 
regions c and d. 



we can obtain A^^^ by solving the fifth-order equation in 
Eq. However, Ap.s can also be obtained by plot- 

ting the left 1l(Ap_s) and right 1r(Ap^s) hand sides of 
Eq. ((25| as functions of Ap^^ , respectively. The real roots 
of Eq. (^5]) can be represented by the crossing points of 
two curves corresponding to Yl(Ap_s) and 1r(Ap^s). Here, 
for clarity, these are shown 



>L(Ap,,) 



G(Ap,,) 



7o« 



(Ap,, - Ao) . 



(29) 
(30) 



The curve and the lines, corresponding to Yl(Ap_s) and 
Yr(Ap^), are schematically shown in Fig. [21 This di- 
agram can be used to analyze the solution of the self- 
consistent equation in Eq. ([^5]) . In Fig.[21Ja), the double- 
well like curve shows how Yl(Ap.s) changes with Ap^^, 
and the straight line shows how Yr(Ap_s) changes with 
Ap^s for a given k, but for larger Aq = Ua — ojq. Fig- 
ure 2(b) shows the solution of the self-consistent equa- 
tion in Eq. (|25p when Aq is in the region inside Fig. 2(a) 
marked by the red dashed box. 

In Fig. [21 the black curves denote the function 
1l(Ap_s)- The functions Yr(Ap_s) are shown by the lines 
with different colors. The slope of the lines in Fig. 2(b) 
are proportional to n. The straight lines in Fig. [2jb) 
correspond to different values of Aq. Recall that Aq 
is the atom-cavity detuning when the two mirrors are 
fixed. Each parameter region for Ag can have at most five 
steady-state solutions and three stable solutions. This 
system has eight parameters. However, during most of 
this study, we will be varying the atom-cavity detun- 
ing Ao (for fixed mirrors), and the elastic constant k. It 
is noted that the different crossing points between the 
transverse axis and the different lines corresponding to 
Yr(Ap s) represent the different detunings Aq. The re- 
gions {a-d) are shown in the figure. The points of in- 
tersection of Fl(Ap,s) and Yr(Ap_s) show the solutions of 
Eq. (|^ for the whole Ap_s region. However, physically, 
it is important to study the solution near the atom-cavity 
detuning Ao, and the unstable solution shown by the yel- 
low dot in Fig. [21(a) can be neglected. 

The letters a,b,c,d at the bottom of Fig. [H^b) repre- 
sent regions with different number of solutions. Below, 
in our discussions for the steady-state solutions, we only 
consider the crossing points between two curves corre- 
sponding to yL(Ap^) and lr{(Ap,s), in the red square 
in Fig. mja). Inside this red square, the number of real 
roots of Eq. can be characterized by three critical 

and A|^'''', when 



values of the detuning Aq: \ ^th ^ 



K]^ < K < oo. When fixing the other parameters, the 
lower bound kl of k can be fixed as shown in the numer- 
ical calculations below. 

As schematically shown in Fig. ^h), there are four 
regions in the space of "roots" divided by the above three 

> 
< 
< 
< 



boundaries: (a) when the atom-cavity detuning Aq 
there is no real root; (b) when A^'^'^ < Aq 
there always exist two roots; (c) when A^"'''' 



A (cd) 
^th ' 

A (cd) 
^th ' 

Ao < 



A 



A, 



(ab) 
th ' 



(be) 
th ' 



there are four real roots; 



th 

(d) when Aq 



there are two real roots. Also at the threshold 

_ A(cd) 



points for Ao — A^J^"'', A['^\ and A[j°^\ the number of 
steady-state roots is one, three, and three, respectively. 
However, in the case when k < kl, the threshold value 



(ab) 
^th 



does not exist, and there are only two threshold 



values, a['^^ and A]."^' , which divide the Ao-parameter 
space into three regions for the roots of Eq. ([^5|) . In 
this case, the number of roots will be explained below for 
given sets of parameters. 



Abe) 
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atom-cavity detuning Aq 

FIG. 3: (Color Online) The steady state solutions Ap'l (i = 
1, 2, 3, 4) versus the atom-cavity detunings Aq (when the two 
mirrors are fixed) for given parameters, e.g., k = lO'^, Ac = 0, 
coa = 10^ 7o = 10-^ 72 = 10-*, g = 10^ n = 2, and ain = 
10. Hereafter, all the quantities are measured in units of 7j^. 
Note that here eight parameters determine the system. Recall 
that Ap,s is the steady-state atom-cavity detuning when one 
mirror is moveable. The inset shows the steady-state solution 
Ap|s, which corresponds to a very large displacement of the 
mirror. As schematically shown in Fig. [5{b), the dotted blue, 
dashed red, and continuous green lines denote the solutions 
Ap^s, Ap^s , and A^p l, respectively. For example, the green 
cross in Fig. 2(b) corresponds to a single value of Aq. Here, Aq 
is swept, and the green cross in Fig. 2(b) becomes a continuous 
curve. Same for one red dot and one blue square in Fig. 2(b). 

We now simulate the four steady state solutions 

aU) A (particular solution label) / ■ rj o A\ 

^p,s = ^probe, steady-state J I* — ^7 4j 

for given parameters as in Ref. [25j, e.g., Ac — 0, uja — 
10^ 7o = 10-^ 72 = 10-^, gVN = 10^, = 2 and 
flin — 10. Hereafter, all quantities are measured in units 
of 7]^. With the above parameters, we can find that the 
lower bound kl is about 6400. 

We first study the case for k < kl = 6400, e.g., k = 

10^. In this case, there are two critical values a['^ = 25 
and A^.'^ = —0.95. Fig. 3 shows the steady-state atom- 
cavity detunings Ap^s (when the mirror moves) versus 
Aq (when the mirror is fixed). We find that there is 
no solution when Aq = Wa — i^o ^ 25. This means that 
when the difference between the atomic frequency cOa and 
the cavity frequency loq is larger than 25, there is no 
steady-state near the detuning Ac. When —0.95 < Aq < 
25, there are two steady-state solutions, i.e., Ap^] and 
A^pJ shown in Fig. [31 In this region, we find that A^pJ 

(2) 

decreases linearly, but Ap^s increases with increasing Aq. 




atom-cavity detuning Aq 

FIG. 4; (Color Online) Steady-state positions (s)!*' {i = 
1, 2, 3, 4) for the moveable mirror versus the atom-cavity de- 
tuning Ao, for the same parameters listed in Fig.O Here, the 
continuous black, dotted blue, dashed red, continuous green 
lines denote (a;)^"'"', {x)''^\ {x)i^\ and {x)i'^\ when the atom- 
cavity detunings Ap,s take the following steady-state values 
Ap^s, A^l, Ap^s, and Ap*s, respectively, as shown in Fig. [S] 

When Aq < —0.95, as shown in Fig. [31 there are four 
steady-state solutions, i.e., Ap'i {i = 1, 2, 3, 4). Figure[31 
shows that when the detuning Aq passes through —0.95, 

/ri\ 

from the right to the left, two additional solutions (ApJ 
and Ap'^]) appear. We also find that two solutions (Ap^s 
and Ap*]) gradually approach Ac = to realize the two- 

(3) 

photon-resonance condition. Moreover, note that Ap^i 
increases almost linearly with increasing Aq. 

With the same parameters as those in Fig. 3, we have 
also plotted in Fig. [4l the Ap-dependent location (a;)^*-* 
of the mirror corresponding to Ap's. We find that there 
is no solution for the steady-state value of (a;)^''' when 
Ao > 25. When -0.95 < Aq < 25, the mirror's position 
(x) , corresponding to the solution Ap^l , exhibits a very 
large (compared with {x)^p , {x)^^\ displacement 
of the mirror, and (a;)^^'' increases as Aq increases. The 
mirror's displacement {x)]. , corresponding to the solu- 

(2) 

tion Ap_s, is nearly zero. When Aq < —0.95, the four 
steady-state solutions for the displacement (x)^ exist si- 
multaneously. Two of them, i.e., (a;)^^-* and {x)f\ show 
that the spring is compressed, and their displacements 
linearly increase when Aq decreases. One of them, i.e., 
{x)^f\ nearly vanishes. 

From Eq. (|16p , we know that the oscillating mirror can 
affect the two-photon resonance by changing the effective 
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frequency lol of the cavity field. Because when the mirror 
is fixed, the two-photon resonant condition becomes 

Ao = C^a - ^0 = Ac. (31) 

However, this condition is modified to 

Ap^s = t^a - = Ac , (32) 

when the mirror is osciUating. We note that the two- 
photon resonant condition in Eq. ()32|) is further modified 
to 

A« = - = Ac , (33) 

when the system reaches its steady state. We assume 
that the two-photon resonant condition in Eq. ([3T|) is 
satisfied when the mirror is fixed. Recall that Ac = 
when Fig. [3] and Fig. |4] are plotted. This means that the 
two-photon resonant condition is Ap,l — in this case. 
We find that {x) ^ is nearly zero as shown in Fig. 2] when 
Aq > Ac = 0. In this case, there is no value of Ap ^ 
approaching zero, as shown in Fig. [3l so the two-photon 
resonance cannot happen. From Fig. [3] and Fig. [H we 
find that the two-photon resonance Ap*s = Ac = might 
happen when Aq < —0.95. Because in this region, ApJ 
and Ap^s can approach zero as in Fig. [31 which correspond 
to the steady-state values of the mirror's positions {x)''p 

and {x)^^\ respectively, as in Fig.lH 

We now turn to study the steady-state values of Ap ^ 
and {x) ^ for the case when 6400 S k < oo, e.g., k — 10'*. 
In this schematically shown in Fig. [Ifb), there 

are three critical values A[|^^ (i — ab, be, cd): ^[if^ = 

2.5 X 10^, aII^'^^ = -0.22 and A^f ^ = -2.5 x 10^ The 
number of solutions for Ap ^ has the same descriptions as 
for Fig.d^b). Similarly to Fig.[3]and Fig.H we plot Fig.E] 

and Fig.[n]for Ap'^ and their corresponding {x)^^^\ respec- 
tively. The two-photon resonance condition in Eq. (j33p 
for Ac = might also be satisfied in this case. Because 

(2) 

when Ao < —0.22, two steady-state values (Ap_s and 

A^S) of Ap^s are near zero, as shown in Fig.[5l which cor- 
respond to the steady-state solution of the mirror's posi- 
tion (a;)^^^ and {x)[^\ as shown in Fig.[6l It is also found 

(2) 

that one, i.e., (x)]. ', of the steady-state solutions of {x) ^ 
nearly vanishes, as shown in Fig. [6l when Ag > Ac = 0. 

Finally, we note that there is only one steady state so- 
lution Ap s = Aq for Eq. ()25|) in the limit k, oo. This 
implies that when the elastic coefficient Mui'j^j is very 
large, it is difficult for the photon pressure to make the 
mirror to even have a tiny displacement, and the oscillat- 
ing mirror does not affect the optomechanical system [30l | . 

IV. ELECTROMAGNETICALLY INDUCED 
TRANSPARENCY IN THE OPTOMECHANICAL 
SYSTEM 

To explore how the mirror's oscillation affects the FIT, 
let us now study the susceptibility of the atomic medium. 




-6-3 3 6 
atom-cavity detuning Aq 

FIG. 5: (Color Online) Steady-state solutions A^^^''^^ {i — 
1, 2, 3, 4) schematically shown in Fig. [2jb) versus Aq with 
K = 10*. The other parameters here are the same as in 
Fig. [3l The inset shows the steady state solution that 
corresponds to a very large displacement of the mirror. As 
shown in Fig. [2jb), the dotted blue, dashed red, and continu- 
ous green lines denote the solutions A^^'^^.j, Ap^|,j , and A^'||,j, 
respectively. 




atom-cavity detuning Aq 

FIG. 6: (Color Online) Steady-positions {x)i'^ {i = 1, 2, 3, 4) 
of the moveable mirror versus the atom-cavity detuning Ao, 
with the same parameters as in Fig. [5] Here, the continu- 
ous black, dotted blue, dashed red, continuous green lines de- 
note {x)i^\ {^)^^\ and {x)i'^\ when the atom-cavity 
detunings Ap(s) (for a moving mirror) takes the following 
steady-state values: A^J^, A^^^^, Ap^s, and A^^l, respectively, 
as shown in Fig. \5\ 



As discussed above, we consider the single-mode cavity 
field as the probe field. In this case, we have the electric 
field 



E{t) = 



H.c. = ee 



— ilul i 



H.c. . (34) 



The linear response of the atomic ensemble to the weak 
probe field E{t) can be described by the susceptibility 



X 



(s) eo ' 



(35) 



Here, the average polarization (p) of the atomic ensemble 
is 



N 



V 



(36) 



i=i 



Using Eqs. (HH), ([M)) and ([55)) . we obtain the suscepti- 
bility X 



X = F 



72S- (Ap,, - Ac)e 



-iF 



S2 + 62 
726 + (Ap,, - Ae) 5 
52 + 92 



with 



F = 



7i (Ap,5 - Ac 



(37) 

(38) 
(39) 



(where F is proportional to the density N/V) and 

e = n^~ Ap,, (Ap,, - AJ + 7i72- (40) 

It is well known that the real, Re(x), and imaginary, 
Ini(x) , parts of the susceptibility x describe the disper- 
sion and absorption of light, respectively. 

As discussed in the last section, one of the solutions 



(1) 



and Ap^l is unstable, and another 



,(4) 



,p_s, is similar to the 
in the region (A^''"' 



^th> There- 



corresponding to x 
solution, corresponding to xi^^ and A 

(2) 

solution Xg and A 
fore, we only consider the other two solutions below. In 
Fig. El Re(xi) and Im(xi) (« = 1, 2) versus Aq are plot- 
ted for the two solutions A^l and Ap'^] studied in Fig. 3. 
Here, Xi and X2 denote the susceptibilities correspond- 



(2) 

ing to Ap/s and A^_'s, respectively. All parameters in 
Fig. 7 are the same as those in Fig. [3l It is found that 
when Aq > Ac = 0, the curves for Rc(xi) and Im(xi ) 
are similar to those in the usual KIT phenomenon [30|. 

This means that the steady state value {x) ]. of the mir- 
ror's displacement is near zero in this region of parame- 
ters, and the oscillating mirror has no effect on the FIT. 

However, when Aq < Ac = 0, the mirror's displace- 
f 2) 

ment {x)\. makes the two-photon-resonance condition 
'21 

U!a — i^L ~ Ac = be approximately satisfied. 



.(3) 
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FIG. 7: (Color Online) The susceptibilities Xi ^i^d X2 
the detuning between the atom and the control field Ac — 
and K = 10^. We plot only two susceptibilities because there 
are three stable solutions in this parameter region but two 
of them are quite similar. Here, the dotted and solid curves 
correspond to the imaginary and real parts, respectively. The 
blue and the green colors correspond to Xi S'Hd the dark red 
and brown colors correspond to X21 respectively. The EIT-like 
region is located between the two peaks of ImXi and Imx2 • 



D 
O 

1/1 

=i 

(/I 

O 



o 



0.02 



0.01 



1 1 

• > 

/ \ 

• i 


1 

1 
1 
1 

1 
1 
1 
1 
1 
■ 

V ■ 
\ 1 
\ 1 

i J 


1 

ilnWi 

1 
1 
1 

■ 
1 
1 

*'•■■■• 


1 1 ^ 


1 


Re^i 



-0.01 

-6-3 3 6 
atom-cavity detuning Aq 

FIG. 8: (Color Online) The susceptibilities Xi and X2 
Ac — 2 and k — 10^. Here, the dotted and solid curves 
correspond to the imaginary and real parts, respectively. The 
blue and the green colors correspond to Xi and the dark red 
and brown colors correspond to X21 respectively. The EIT-like 
region is located between the two peaks of ImXi and Imx2- 
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FIG. 9; (Color Online) The susceptibilities for Ac = and 
K — 1.6 X 10^". Note that this value of the elastic constant 
K is huge, corresponding to an almost fixed mirror. In this 
case, the ReXi and Imxi occur for mostly positive values 
of the atom-cavity detuning. When the spring constant be- 
comes softer, as in Fig. 7, the ReXi and Imxj have a response 
that extend over a huge range of values of the atom-cavity 
detuning Aq, even for Ac > —6. When k tends to infinity, 
the susceptibility becomes the same as that in the usual EIT 
phenomenon. This consistency check is reassuring for our 
calculations. Here, the dotted and solid curves correspond to 
the imaginary and real parts, respectively. The blue and the 
green colors correspond to X\ a^nd the dark red and brown col- 
ors correspond to Xi i respectively. The EIT-like region occurs 
near zero detuning Aq. 



As a result, in this region, each of the curves Re(xi) 
and Im(x]^) is almost close to zero, like an infinite "tail". 
When Aq < —0.95 as shown in Fig. [31 another solution 



f 3 o 

ApJ emerges, thus we also have a Ap'^i-dependent sus- 
ceptibihty X2- 

As shown in Fig. [3 the curves corresponding to the real 
Re(x2) E^nd imaginary Im(x2) parts of the susceptibility 
X2 are similar to those in the usual EIT [s^l, since the 



(3) 



mirror's displacement 



,(3) 



is nearly zero. From Fig. [3 



we find that the right part of the curve Im(x2) almost 
merges with the left part of the curve Imxi, so that a 
transparency "window" is formed. 

To know how Ac affects the EIT, we can also study 
the susceptibility x for the detuning Ac = 2 when other 
parameters are assumed to be the same as those in Fig. [T] 
In this case, the steady-state solutions of Ap^s and {x)^ 
are similar to those for Ac = 0. The curves corresponding 
to these solutions are just rightward shifts for the curves 
in Fig. [3] and Fig. [H but the shapes of the curves are al- 
most the same. Similar to Fig. [3 we choose two steady- 
state solutions and plot Re(xi) and Im(xi) {i = 1, 2) 



versus Aq. Obviously, the essential conclusions remain 
unchanged as those in Fig. [3 but all curves have a right- 
ward shift. 

Based on the analysis in Sec. Ill, a[|^^'' and a|JJ'''' ap- 
proach zero from the left side when k is increased. When 
K > kl, e.g., K > 6400 in Fig. [5] and Fig.[6l the larger k 
corresponds to shorter "tails" of the curves correspond- 
ing to Re(xi) and Im(xi) when Aq > —0.95. In the 
limit oo, the "tails" (for Aq > 0) disappear and the 
right part of the curve Im(x2) just meets the left part 
of the curve Im(xi) to form a transparency "window". 
In this limit, all the physical properties return to that 
in the usual EIT phenomenon. Namely, we recover the 
standard EIT when k — > oo. This asymptotic result is 
shown in Fig. [5] with a large k (e.g., k = 1.6 x lO"'^"), but 
other parameters are the same as in Fig. [31 In Fig. [51 
around the point Aq > Ac = 0, the left parts of the 
curves Re(xi) and Im(xi) nearly merge with the right 
parts of the curves Re(x2) and Im(x2)7 and thus trans- 
parency windows are formed as in the usual EIT [30| . 



V. CONCLUSION AND REMARKS 

We have studied the effects of the end mirror's oscil- 
lation on the EIT phenomenon for an atomic ensemble 
confined in a gas cell placed in a micro-cavity. This study 
can help us to quantitatively consider the quantum inter- 
face between an optomechanical system and an atomic 
gas displaying EIT. The results obtained could be used 
to improve the measurement precision based on the EIT 
effect, when the medium is placed inside a microcavity 
with a one-end oscillating mirror and the cavity acts as 
the probe light. We have shown that the whole system 
exhibits multi-stability due to the mirror's vibration, and 
we have also numerically obtained the threshold values 
of the parameters, which can help determine how many 
steady-states exist in the corresponding parameter re- 
gions. This multi-stability can be explicitly displayed 
through a modified EIT phenomenon. Consequently, we 
investigate the effects of the multi-steady-state solutions 
on the EIT phenomenon and find that in some param- 
eter regions there are two solutions that approximately 
satisfy the two-photon resonance condition. Therefore, 
the properties of both the real and imaginary parts of 
the susceptibility are significantly altered. When the 
spring elastic constant k increases, the mirror becomes 
less moveable, and in this case all properties of the system 
gradually revert to those of the usual EIT phenomenon. 
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